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Abstract: We revisit the calculation of the relic density of the lightest Kaluza-Klein 
particle (LKP) in the model of Universal Extra Dimensions. The Kaluza-Klein (KK) 
particle spectrum at level one is rather degenerate, and various coannihilation processes 
may be relevant. We extend the calculation of hep-ph/0206071 to include coannihilation 
processes with all level one KK particles. In our computation we consider a most general 
KK particle spectrum, without any simplifying assumptions. In particular, we do not 
assume a completely degenerate KK spectrum and instead retain the dependence on each 
individual KK mass. As an application of our results, we calculate the Kaluza-Klein relic 
density in the Minimal UED model, turning on coannihilations with all level one KK 
particles. We then go beyond the minimal model and discuss the size of the coannihilation 
effects separately for each class of level 1 KK particles. Our results provide the basis for 
consistent relic density computations in arbitrarily general models with Universal Extra 
Dimenions. 
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1. Introduction 

The Standard Model (SM) of particle physics has been astonishingly successful in explaining 
much of the presently available experimental data. However, it still leaves open a number of 
outstanding fundamental questions whose answers are expected to emerge in a more general 
theoretical framework which extends the SM at higher energy scales. The main motivation 
for pushing the energy frontier beyond the Terascale comes from two main issues. The first 
is the question, what breaks the electroweak gauge symmetry, and if it is a Higgs field, 
why is the Higgs particle so light. The second is the dark matter problem, which has no 
explanation within the SM. By now we have accumulated a wealth of astrophysical data 
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which all point to the existence of a dark, non-baryonic component of the matter in the 
Universe. The most recent WMAP data confirm the standard cosmological model and pin 
down the amount of cold dark matter 1 as 0.094 < QcDMh 2 < 0.129, which is consistent 
with earlier indications, but much more precise. However, the microscopic nature of dark 
matter is presently unknown, as all known particles are excluded as dark matter candidates. 
This makes the dark matter problem the most pressing phenomenological motivation for 
particles and interactions beyond the Standard Model. 

Among the most attractive explanations of the dark matter problem is the "WIMP" hy- 
pothesis - that dark matter consists of stable, neutral, weakly interacting massive particles. 
Indeed, many theoretically motivated extensions of the Standard Model which attempt to 
resolve the gauge hierarchy problem in the electroweak sector, already contain particles 
which can be identified as WIMP dark matter candidates. A straightforward calculation 
of the WIMP thermal relic abundance, which we shall review in some detail in Section |3|, 
reveals that WIMP particles with masses in the TeV range may provide most of the dark 
matter. 

Supersymmetry (SUSY) and Extra Dimensions, which appear rather naturally in string 
theory, are the leading candidate theories for physics beyond the SM. By now they are both 
known to possess WIMP dark matter candidates. The collider and astroparticle signals 
of SUSY dark matter have been extensively studied [1,2]. Recently, WIMPs from extra 
dimensional theories, have started to attract similar interest. A particularly appealing sce- 
nario is offered by the so called Universal Extra Dimensions (UED), originally proposed 
in [3], where all SM particles are allowed to freely propagate into the bulk of one or more 
extra dimensions. The case of UED bears interesting analogies to supersymmetry, and 
sometimes has been referred to as "bosonic supersymmetry" [4]. Many of the virtues of 
supersymmetry remain valid in the UED framework. For example, the existence of a dark 
matter candidate, the lightest Kaluza-Klein particle (LKP), is guaranteed by a discrete 
symmetry, called KK-parity. KK-parity also eliminates dangerous tree-level contributions 
to electroweak precision observables, rendering the model viable for a large range of pa- 
rameters below the TeV scale. 

There are also some important differences between SUSY and UED. For example, 
the mass spectrum of Kaluza-Klein particles is rather degenerate, even after radiative 
corrections [5]. This means that the computation of the LKP relic density gets rather 
complicated, since there are many possible coannihilation processes with other particles in 
the KK spectrum [6]. This situation should be contrasted with the case of supersymmetric 
models where typically the spectrum is widely spread 2 , and barring fine-tuned coincidences, 
only one or at most a few additional particles participate in coannihilation processes. 

A second important distinction between SUSY and UED is the presence of a KK tower 
in the case of extra dimensions. The KK masses roughly scale as n/R, where R is the size 

1 From now on and throughout the paper we shall use Q = Qcdm to denote the dark matter relic density. 

2 Notice that the recently proposed little Higgs models with T-parity [7-9] are reminiscent of UED, and 
their spectrum does not have to be degenerate, so they are four dimensional examples which bear analogies 
to both UED and SUSY. 
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of the extra dimension, and n is the KK level. This has an important implication for the 
KK relic density calculation, since the n = 2 particles naturally have masses twice as large 
as the LKP mass, and resonant annihilation processes may become relevant [10,11]. 

A third important difference is encoded in the spins of the new particles. The su- 
perpartners have spins which differ by 1/2 unit from their SM counterparts, while the 
spins of the KK particles are the same as in the SM. This has important consequences for 
the dark matter abundance as well. For example, the dark matter candidate in SUSY is 
usually the lightest neutralino, which is in general some mixture of a Bino, a Wino and 
Higgsinos (the superpartners of the hypercharge gauge boson, the neutral SU(2)w gauge 
boson, and the neutral Higgs bosons, respectively). The neutralino is a Majorana particle, 
and its annihilation to SM fermion final states is helicity suppressed. In contrast, the dark 
matter candidate in UED is the KK mode of a Higgs or gauge boson, and its annihilation 
cross-sections have no anomalous suppressions. For example, in the Minimal UED model 
(MUED), which we review in more detail below in Section ^, the LKP is B\, the KK part- 
ner of the hypercharge gauge boson, and is a spin one particle. The preferred mass range 
for the LKP is therefore somewhat larger than in super symmetry. 

The first and only comprehensive calculation of the UED relic density to date was per- 
formed in [6]. The authors considered two cases of LKP: the KK hypercharge gauge boson 
B\ and the KK neutrino v\. The case of B\ LKP is naturally obtained in MUED, where the 
radiative corrections to B\ are the smallest in size, since they are only due to hypercharge 
interactions. The authors of [6] also realized the importance of coannihilation processes 
and included in their analysis coannihilations with the S'[/(2)vi/-singlet KK leptons, which 
in MUED are the lightest among the remaining n = 1 KK particles. It was therefore 
expected that their coannihilations will be most important. Subsequently, Refs. [10, 11] 
analyzed the resonant enhancement of the n = 1 (co) annihilation cross-sections due to 
n = 2 KK particles. 

Our goal in this paper will be to complete the LKP relic density calculation of Ref. [6]. 
We will attempt to improve in three different aspects: 

• We will include coannihilation effects with all n = 1 KK particles. The motivation 
for such a tour de force is twofold. First, recall that the importance of coannihilations 
is mostly determined by the degeneracy of the corresponding particle with the dark 
matter candidate. In the Minimal UED model, the KK mass splittings are due almost 
entirely to radiative corrections. In MUED, therefore, one might expect that, since 
the corrections to KK particles other than the KK leptons are relatively large, their 
coannihilations can be safely neglected. However, the Minimal UED model makes 
an ansatz [5] about the cut-off scale values of the so called boundary terms, which 
are not fixed by known SM physics, and are in principle arbitrary. In this sense, the 
UED scenario should be considered as a low energy effective theory with a multitude 
of parameters, just like the MSSM, and the MUED model should be treated as 
nothing more than a simple toy model with a limited number of parameters, just like 
the "minimal supergravity" version of supersymmetry, for example. If one makes a 
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different assumption about the inputs at the cut-off scale, both the KK spectrum 
and its phenomenology can be modified significantly. In particular, one could then 
easily find regions of this more general parameter space where other coannihilation 
processes become active. On the other hand, even if we choose to restrict ourselves 
to MUED, there is still a good reason to consider the coannihilation processes which 
were omitted in the analysis of [6] . While it is true that those coannihilations are more 
Boltzmann suppressed, their cross-sections will be larger, since they are mediated by 
weak and/or strong interactions. Without an explicit calculation, it is impossible to 
estimate the size of the net effect, and whether it is indeed negligible compared to 
the purely hypercharge-mediated processes which have already been considered. 

• We will keep the exact value of each KK mass in our formulas for all annihilation 
cross-sections. This will render our analysis self-consistent. All calculations of the 
LKP relic density available so far [6,10,11], have computed the annihilation cross- 
sections in the limit when all level 1 KK masses are the same. This approximation 
is somewhat contradictory in the sense that all KK masses at level one are taken to 
be degenerate with LKP, yet only a limited number of coannihilation processes were 
considered. In reality, a completely degenerate spectrum would require the inclusion 
of all possible coannihilations. Conversely, if some coannihilation processes are being 
neglected, this is presumably because the masses of the corresponding KK particles 
are not degenerate with the LKP, and are Boltzmann suppressed. However, the 
masses of these particles may still enter the formulas for the relevant coannihilation 
cross-sections, and using approximate values for those masses would lead to a certain 
error in the final answer. Since we are keeping the exact mass dependence in the 
formulas, within our approach heavy particles naturally decouple, coannihilations 
are properly weighted, and all relevant coannihilation cross-sections behave properly. 
Notice that the assumption of exact mass degeneracy overestimates the corresponding 
cross-sections and therefore underestimates the relic density. This expectation will 
be confirmed in our numerical analysis in Section ||. 

• We will try to improve the numerical accuracy of the analysis by taking into account 
some minor corrections which were neglected or approximated in [6]. For example, we 
will use a temperature-dependent g* (the total number of effectively massless degrees 
of freedom, given by eq. fl3.6|) below) and include subleading corrections fl3.19|) in the 
velocity expansion of the annihilation cross-sections. 

The availability of the calculation of the remaining coannihilation processes is impor- 
tant also for the following reason. Coannihilations with SU {2)w-sing\et KK leptons were 
found to reduce the effective annihilation cross-section, and therefore increase the LKP 
relic density. This has the effect of lowering the range of cosmologically preferred values 
of the LKP mass, or equivalently, the scale of the extra dimension. However, one could 
expect that coannihilations with the other n = 1 KK particles would have the opposite 
effect, since they have stronger interactions compared to the S , C/(2)vy-singlet KK leptons 
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and the B\ LKP. As a result, the preferred LKP mass range could be pushed back up. 
For both collider and astroparticle searches for dark matter, a crucial question is whether 
there is an upper limit on the WIMP mass which could guarantee discovery, and if so, 
what is its precise numerical value. To this end, one needs to consider the effect of all 
coannihilation processes which have the potential to enhance the LKP annihilations. We 
will see that the lowering of the preferred LKP mass range in the case of coannihilations 
with S'C/(2)vy-singlet KK leptons is more of an exception rather than the rule, and the 
inclusion of all remaining processes is needed in order to derive an absolute upper bound 
on the LKP mass. 

The paper is organized as follows. In Section || we review the Minimal model of Uni- 
versal Extra Dimensions (MUED). In Section |3.l| fl3.2|) we review the standard calculation 
of the WIMP relic density in the absence (presence) of coannihilations. Then in Section [| 
we present our results for the LKP relic abundance in MUED. In Section || we extend 
our analysis beyond MUED and investigate the size of the coannihilation effects for each 
class of KK particles at level 1. We summarize and conclude in Section |6|. The Appendix 
contains a list of our formulas for all relevant annihilation cross-sections, in the limit of 
equal masses. Those are given only for reference and contact with previous work [6, 12], 
since in our numerical code we use the corresponding different values for the individual 
masses. 

2. Universal Extra Dimensions 

The simplest version of UED has all the SM particles propagating in a single extra di- 
mension of size R, which is compactified on an S1/Z2 orbifold. More complicated models 
have also been considered, motivated by ideas about electroweak symmetry breaking and 
vacuum stability [13-15], neutrino masses [16,17], proton stability [18], the number of gen- 
erations [19] or fermion chirality [20,21]. A peculiar feature of UED is the conservation of 
Kaluza-Klein number n at tree level, which is a simple consequence of momentum conser- 
vation along the extra dimension. However, bulk and brane radiative effects [5,22,23] break 
KK number down to a discrete conserved quantity, the so called KK parity, (— l) n . KK 
parity ensures that the lightest KK partners - those at level one - are always pair-produced 
in collider experiments, similar to the case of supersymmetry models with conserved im- 
parity. On the one hand, this leads to rather weak bounds on KK partner masses from 
direct searches at colliders [3,4]. On the other hand, the collider signatures of UED are 
very similar to those of supersymmetry, and discriminating between the two scenarios at 
lepton [24-28] or hadron [29-36] colliders is currently an active field of study. KK parity 
conservation also implies that the contributions to various precisely measured low-energy 
observables [37-47] only arise at loop level and are small. As a result, the limits on the 
scale of the extra dimension, from precision electro-weak data, are rather weak, constrain- 
ing R^ 1 to be larger than approximately 250 GeV [41]. An attractive feature of UED is the 
presence of a stable massive particle which can be a cold dark matter candidate [5,6,48,50]. 
The lightest KK partner (LKP) at level one is also the lightest particle with negative KK 
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parity and is stable on cosmological scales. The identity of the LKP is a delicate issue, 
however, as it depends on the interplay between the one-loop radiative corrections to the 
KK mass spectrum and the brane terms generated by unknown physics at high scales [5]. 
In the Minimal UED model defined in [4], the LKP turns out to be the KK partner B\ 
of the hypercharge gauge boson [5] and its relic density is typical of a WIMP candidate: 
in order to explain all of the dark matter, the Bi mass should be in the range 600 — 800 
GeV, depending on the rest of the KK spectrum [6,10,11,49]. The experimental signals of 
Kaluza-Klein dark matter have also been discussed and it has been realized that it offers 
excellent prospects for direct [50-52] or indirect detection [50,53-61]. 

In the Minimal UED model, the bulk interactions of the KK modes readily follow from 
the SM Lagrangian and contain no unknown parameters other than the mass, rrih, of the 
SM Higgs boson. In contrast, the boundary interactions, which are localized on the orbifold 
fixed points, are in principle arbitrary, and their coefficients represent new free parameters 
in the theory. Since the boundary terms are renormalized by bulk interactions, they are 
scale dependent [22] and cannot be completely ignored since they will be generated by 
renormalization effects. Therefore, we need an ansatz for their values at a particular scale. 
As with any higher dimensional Kaluza-Klein theory, the UED model should be treated 
only as an effective theory which is valid up to some high scale A, at which it matches to 
some more fundamental theory. The Minimal UED model therefore has only two input 
parameters: the size of the extra dimension, R, and the cutoff scale, A. The number of 
KK levels present in the effective theory is simply AR and may vary between a few and 
~ 40, where the upper limit comes from the breakdown of perturbativity already below 
the scale A. Unless specified otherwise, for our numerical results below, we shall always 
choose the value of A so that AR = 20, although analyses of unitarity constraints in the 
higher dimensional Standard Model typically yield a lower bound on A [62,63]. Changing 
the value of A will have very little impact on our results since the A dependence of the KK 
mass spectrum is only logarithmic. For the SM Higgs mass rrih we shall adopt the value 
m h = 120 GeV. 

3. The Basic Calculation of the Relic Density 
3.1 The standard case 

We first summarize the standard calculation for the relic abundance of a particle species 
X which was in thermal equilibrium in the early universe and decoupled when it became 
nonrelativistic [6,64,65]. The relic abundance is found by solving the Boltzmann equation 
for the evolution of the \ number density n 

_ = _ 3ffn _ (cJ7; )( n 2_ n 2 g) (31) 

where H is the Hubble parameter, v is the relative velocity between two %'s, (<rv) is the 
thermally averaged total annihilation cross-section times relative velocity, and n eq is the 
equilibrium number density. At high temperature (T S> m), n eq ~ T 3 (there are roughly 
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as many x particles as photons). At low temperature (T <C m), in the nonrelativistic 
approximation, n eq can be written as 



.'/(^r)'' 7 • (3.2) 



/ mT 



where m is the mass of the relic x> T 1S the temperature and g is the number of internal 
degrees of freedom of x such as spin, color and so on. We see from eq. (|3.2[) that the density 
n eq is Boltzmann-suppressed. At high temperature, x particles are abundant and rapidly 
convert to lighter particles and vice versa. But shortly after the temperature T drops below 
m, the number density decreases exponentially and the annihilation rate T = (av)n drops 
below the expansion rate H. At this point, x' s stop annihilating and escape out of the 
equilibrium and become thermal relics, (av) is often approximated by the nonrelativistic 



expansion 3 



where 



(av) =a + b(v 2 ) + 0({v 4 )) « a + 6b/x + O { — ) , (3.3) 



m 

x = - . (3.4) 

By solving the Boltzmann equation analytically with appropriate approximations [6,64,65], 
the abundance of x is given by 

^ l2 1.04 x 10 9 x F 1 . . 

0,/i 2 » — , (3.5) 

where the Planck mass Mp\ = 1.22 x 10 19 GeV and g* is the total number of effectively 
massless degrees of freedom, 

7 

i=bosons i=fermions 

The freeze-out temperature, x F , is found iteratively from 

i ft , Q^ f^> 9 mM Pl (a + 6b/x F ) \ 
x F = ln c(c + 2)W— — g , 3.7) 

V V 8 27r V9*(x F )x F J 

where the constant c is determined empirically by comparing to numerical solutions of the 
Boltzmann equation and here we take c = ^ as usual. The coefficient | in the right hand 
side of (|3.6| ) accounts for the difference in Fermi and Bose statistics. Notice that g* is 
a function of the temperature T, as the thermal bath quickly gets depleted of the heavy 
species with masses larger than T. 



3 Note, however, that the method fails near s-channel resonances and thresholds for new final states [66]. 
In the interesting parameter region of UED, we are always sufficiently far from thresholds, while for the 
treatment of resonances, see [10,11]. 
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3.2 The case with coannihilations 



When the relic particle % is nearly degenerate with other particles in the spectrum, its 
relic abundance is determined not only by its own self-annihilation cross-section, but also 
by annihilation processes involving the heavier particles. The previous calculation can be 
generalized to this "coannihilation" case in a straightforward way [6,65,66]. Assume that 
the particles Xi are labeled according to their masses, so that m,j < mj when i < j. The 
number densities rij of the various species Xi obey a set of Boltzmann equations. It can be 
shown that under reasonable assumptions [66] , the ultimate relic density n of the lightest 
species xi (after all heavier particles Xi have decayed into it) obeys the following simple 
Boltzmann equation 

dn 



= -3Hn - (a eff v)(n 2 - n 2 eq ) , (3.8) 



dt 



where 



a eff (x) = X;^S^(l + A i ) 3 / 2 (l + A j ) 3 / 2 exp(-x(A i + A J )) , (3.9) 



N 



9eff 



9eff(x) = 9i0- + Ai) 3/2 exp(-xA l ) , (3.10) 
i=i 

Aj = m ~ mi . (3.11) 
mi 

Here aij = cr(xiXj ~^ SM), gi is the number of internal degrees of freedom of particle Xi 
and n = YliLi n « * s the density of xi we want to calculate. This Boltzmann equation can 
be solved in a similar way [6,66], resulting in 

u2 _ 1.04 x 10 9 x F 1 

~ 7T — — , (3.12) 

m pi ^/gZMla + ZIb/xp 



with 



/•oo 

I a = xp I a e ff(x)x~ 2 dx , (3.13) 

J Xp 



I b = 2x 2 F / b eff (x)x- 3 dx . (3.14) 



CO 



Xp 



The corresponding formula for xf becomes 

XF = In (r(r + 2)^^ mM»(a e// ( x F ) + 66^ ^)/^ 

V » 8 2f y/g*(x F )x F J 

Here a e // and 6 e // are the first two terms in the velocity expansion of <j e // 

<7 e //(x) V = O e //(x) + 6 e //(x) V 2 + 0(v 4 ) . (3.16) 
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Comparing eqs. ( |3.9| ) and ( |3.16| ), one gets 

TV 

«e//(*) = + A0 3/2 (l + A,0 3/2 exp(-x(A i + A,)) , (3.17) 



9eff 



b eff(x) = ^^•^(l + A J ) 3/2 (l + A J ) 3 / 2 exp(- a; (A i + A,)) , (3.18) 

ij 9 eff 

where aij and b%j are obtained from aijV = dij + 6r,t> 2 + 0(v A ). 

Considering relativistic corrections [64, 67, 68] to the above treatment results in an 
additional subleading term which can be accounted for by the simple replacement 

b^b- -a (3.19) 
4 

in the above formulas. 

4. Relic Density in Minimal UED 

For the purposes of our study we have implemented the relevant features of the Minimal 
UED model in the CompHEP event generator [69]. We incorporated all n = 1 and n = 2 KK 
modes as new particles, with the proper interactions and one-loop corrected masses [5]. 
Similar to the SM case, the neutral gauge bosons at level 1, Z\ and 71, are mixtures of 
the KK modes of the hypercharge gauge boson and the neutral SU(2)\y gauge boson. 
However, as shown in [5], the radiatively corrected Weinberg angle at level 1 and higher 
is very small. For example, 71, which is the LKP in the minimal UED model, is mostly 
the KK mode of the hypercharge gauge boson. Therefore, for simplicity, in the code we 
neglected neutral gauge boson mixing for n = 1. We then use our UED implementation 
in CompHEP to derive analytic expressions for the (co) annihilation cross-sections between 
any pair of n = 1 KK particles. Our code has been subjected to numerous tests and cross- 
checks. For example, we reproduced all results from [6]. We have also used the same code 
for independent studies of the collider and astroparticle signatures of UED [25,33,61,70] 
and thus have tested it from a different angle as well. 

The mass spectrum of the n = 1 KK partners in Minimal UED can be found, for 
example, in Fig. 1 of [4]. In MUED the next-to-lightest KK particles are the singlet KK 
leptons and their fractional mass difference from the LKP is 4 

mp„, — 

A tm = — ^ ^ ~ 0.01 . 4.1) 

Notice that the Boltzmann suppression 

e -A tm x F ^ g-0.01-25 = e -0.25 



4 In this paper we follow the notation of [6] where the two types of n = 1 Dirac fermions are distinguished 
by an index corresponding to the chirality of their zero mode partner. For example, Irx stands for an 
S(/(2)iv-singlet Dirac fermion, which has in principle both a left-handed and a right-handed component. 
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is not very effective and co-annihilation processes with Irx are definitely important, hence 
they were considered in [6]. What about the other, heavier particles in the n = 1 KK 
spectrum in MUED? Since their mass splittings from the LKP 



Aj = m (4.2) 

are larger, their annihilations suffer from a larger Boltzmann suppression. However, the 
couplings of all n = 1 KK partners other than £ri are larger compared to those of 71 and 
£ri. For example, S'C/(2)vK-doublet KK leptons £li couple weakly, and the KK quarks 
qi and KK gluon g\ have strong couplings. Therefore, their corresponding annihilation 
cross-sections are expected to be larger than the cross-section of the main 7171 channel. 

We see that for the other KK particles, there is a competition between the increased 
cross-sections and the larger Boltzmann suppression. An explicit calculation is therefore 
needed in order to evaluate the net effect of these two factors, and judge the importance 
of the coannihilation processes which have been neglected so far. One might expect that 
coannihilations with 5C/(2)vi/-doublet KK leptons might be numerically significant, since 
their mass splitting in MUED is ~ 3% and the corresponding Boltzmann suppression factor 
is only e -o.03-25 „ e -o.75. 

In our code we keep all KK masses different while we neglect all the masses of the 
Standard Model particles. As an illustration, let us show the a and b terms for 7171 
annihilation only. For fermion final states we find the a-term and 6-term of £7(7171 — > ff)v 
as follows 



32^^ / Y f L Yf R . 



^4*«lN c mi ( vt 11m,, + 14m m f[i - 13m f 



27 \ jL (m\ + m] 



Um* + Um 2 m 2 f - 13m 4 f 



where 51 is the gauge coupling of the hypercharge (7 (l)y gauge group, a± = |^ and A^ c 
for f = q and A^ c = 1 for / = £ is the hypercharge of the fermion /. 
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Figure 1: The a-term of the annihilation cross-section for (a) 7171 — ► e + e~ and (b) 7171 — ► <fi(j>* , as 
a function of the mass of the t-channel particle(s). We fix the LKP mass at m 11 — 500 GeV and vary 
(a) the KK lepton mass m eR1 — m eL1 or (b) the KK Higgs boson mass m^ >1 . The blue solid lines 



are the exact results ( |4.j|j and ( [4.7| ), while the red dotted lines correspond to the approximations 



(4.4) and (4.8). 



For the Higgs boson final states we get 

" Y 108 <V «+<) 4 ) { ' 

1 = ±± . (4.10) 

In the limit where all KK masses are the same (the second line in each formula above), 
we recover the result of [6]. Notice the tremendous simplification which arises as a result 
of the mass degeneracy assumption. In Fig. |l] we show the a terms of the annihilation 
cross-section for two processes: (a) 7171 — > e + e~ and (b) 7171 — ► 00*, as a function of 
the mass of the t-channel particle(s). We fix the LKP mass at m 7l = 500 GeV and vary 
(a) the KK lepton mass m em = m eL1 or (b) the KK Higgs boson mass m^ x . The blue 
solid lines are the exact results ( |4.3| ) and (|4.7|) , while the red dotted lines correspond to 
the approximations ( |4.4j ) and ( |4.8| ) in which the mass difference between the t-channel 
particles and the LKP has been neglected. We see that the approximations fl4.4|) and (fOI) 
can result in a relatively large error, whose size depends on the actual mass splitting of the 
KK particles. This is why in our code we keep all individual mass dependencies. 
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Figure 2: The quantities (a) g*(Tp) and (b) xp as a function of R 1 in MUED. 

Another difference between our analysis and that of Ref. [6] is that here we shall use 
a temperature-dependent g* function as defined in ([0]). The relevant value of which 
enters the answer for the LKP relic density ( [3.12 ) is g*(Tp), where Tp = m^/xp is the 



freeze-out temperature. In Fig. [2]a we show a plot of g*(Tp) as a function of R 1 in 
MUED, while in Fig. 0b we show the corresponding values of xp. In Fig. |a one can 
clearly see the jumps in g* when crossing the bb, W + W~, ZZ and hh thresholds (from left 
to right). The ti threshold is further to the right, outside the plotted range. As we shall see 
below, cosmologically interesting values of Qh 2 are obtained for below 1 TeV, where 
g*(Tp) = 86.25, since we are below the W + W~ threshold. The analysis of Ref. [6] assumed 
a constant value of <?* =92, which is only valid between the IU + VF _ and ZZ thresholds. 

The expert reader has probably noticed from Fig. ||b that the values of xp which we 
obtain in MUED are somewhat larger than the x f values one would have in typical SUSY 
models. This is due to the effect of coannihilations, which increase g e ff (see Fig. [|c below) 
and therefore xp, in accordance with ( |3.15| ). 

We are now in a position to discuss our main result in MUED. In Fig. || we show 
the LKP relic density as a function of i? -1 in the Minimal UED model. We show the 
results from several analyses, each under different assumptions, in order to illustrate the 
effect of each assumption. We first show several calculations for the academic case of no 
coannihilations. The three solid lines in Fig. |3| account only for the 7171 process. The 
(red) line marked "a" recreates the analysis of Ref. [6], assuming a degenerate KK mass 
spectrum. The (blue) line marked "b" repeats the same analysis, but uses T-dependent g* 
according to ( |3,6j ) and includes the relativistic correction to the 6-term ( 3.1 9| ) . The (black) 



line marked "c" further relaxes the assumption of KK mass degeneracy, and uses the actual 
MUED mass spectrum. 

Comparing lines "a" and "b", we see that, as already anticipated from Fig. |2|a, ac- 
counting for the T dependence in g* has the effect of lowering g^(xp), a e ff(xp), and 
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Figure 3: Relic density of the LKP as a function of R^ 1 in the Minimal UED model. The (red) line 
marked "a" is the result from considering 7171 annihilation only, following the analysis of Ref. [6], 
assuming a degenerate KK mass spectrum. The (blue) line marked "b" repeats the same analysis, 
but uses T -dependent <?* according to (3.6) and includes the relativistic correction to the b-term 
( 3.15p . The (black) line marked "c" relaxes the assumption of KK mass degeneracy and uses the 
actual MUED mass spectrum. The dotted line is the result from the full calculation in MUED, 
including all coannihilation processes, with the proper choice of masses. The green horizontal band 
denotes the preferred WMAP region for the relic density 0.094 < VLcDMh 2 < 0.129. The cyan 
vertical band delineates values of R^ 1 disfavored by precision data. 



correspondingly, increasing the prediction for Qh 2 . This, in turns, lowers the preferred 
mass range for 71. Next, comparing lines "b" and "c", we see that dropping the mass 
degeneracy assumption has a similar effect on a e ff(xp) (see Fig. |]), an d further increases 
the calculated J7/i 2 . This can be easily understood from the t-channel mass dependence 
exhibited in ( ^4.3[ ) and ( |4.7| ). The t-channel masses appear in the denominator, and they 
are by definition larger than the LKP mass. Therefore, using their actual values can only 
decrease c e // and increase Qh 2 . 

The dotted line in Fig. ^ is the result from the full calculation in MUED, including 
all coannihilation processes, with the proper choice of masses. The green horizontal band 
denotes the preferred WMAP region for the relic density 0.094 < VL C DNih 2 < 0.129. The 
cyan vertical band delineates values of R^ 1 disfavored by precision data [41]. We see that 
according to the full calculation, the cosmologically ideal mass range is m 7l ~ 500 — 600 
GeV, when 71 accounts for all of the dark matter in the Universe. This range is somewhat 
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lower than earlier studies have indicated, mostly due to the effects discussed above. Since 
the MUED model will be our reference point for the investigations in Section ||, the dotted 
line from Fig. ^ will be appearing in all subsequent plots in Section [| below. 



5. Relative Weight of Different Coannihilation Processes 

As we already explained in the Introduction, the assumptions behind the MUED model 
can be easily relaxed by allowing nonvanishing boundary terms at the scale A. This would 
modify the KK spectrum and correspondingly change our prediction for the KK relic 
density from the previous section. Our code is able to handle such more general cases with 
ease, since we use as inputs the physical KK masses. In order to gain some insight into 
the cosmology of such non-minimal scenarios, we have studied the effects of varying the 
n = 1 KK masses one at a time. The change in any given KK mass will not only enhance 
or suppress the related coannihilation processes, but also impact any other cross-sections 
which happen to have a dependence on the mass parameter being varied. Thus the results 
in this section may allow one to judge the importance of each individual coannihilation 
process, and anticipate the answer for Qh? in non-minimal models. 



We have classified the discussion in this section by particle types. Section ^lj contains 
our results for the annihilation processes with KK leptons. Many of our results have already 
appeared in Ref. [6]. The new element here is the discussion of Eli coannihilations. The 



results presented in Sections 5.2 and 5.3 are completely new - there we investigate the 
coannihilation effects with strongly interacting KK modes and electroweak gauge bosons 
and/or Higgs bosons, respectively. 

5.1 Effects due to coannihilations with KK leptons 

We begin with a discussion of 71 coannihilations with the n = 1 S77(2) ^-singlet leptons 
Irx and the n = 1 SU(2) ^-doublet leptons £li- One might expect that those processes 
will be important, since the KK leptons receive relatively small one-loop mass corrections. 
For example, in the Minimal UED model Ae R1 ~ 1% an d &£ L1 ~ 3%. It is natural to 
expect that this degeneracy might persist in non-minimal models as well. 

Our approach is as follows. Since we keep separate values for the KK masses, when 
we start varying any one of them, we have to somehow fix the remainder of the KK mass 
spectrum. We choose to use MUED as our reference model, hence the masses which are not 
being varied, will be fixed according to their MUED values. We shall still show results for 
£lh 2 as a function of R , but for various fixed values of the corresponding mass splitting 
Aj defined in eq. (|4.2j ). We shall also always display the reference MUED model line, for 
which, of course, Aj takes its MUED value. 

Our first example is shown in Fig. ^j, where we illustrate the size of the coannihilation 
effects for (a) 1 generation or (b) 3 generations of degenerate singlet KK leptons £ri. The 
lines show the LKP relic density as a function of R , for several choices of the mass 
splitting Ai R1 between the LKP and the S'C/(2)vi/-singlet KK fermions £r\. The solid lines 
from top to bottom in both (a) and (b) correspond to A ej?1 = 0,0.3,0.1, and the dotted 
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Figure 4: Coannihilation effects of (a) 1 generation or (b) 3 generations of singlet KK leptons. 
The lines show the LKP relic density as a function of R~ l , for several choices of the mass splitting 
A^ H1 between the LKP and the SV "(2) w -singlet KK fermions In each case we use the MUED 
spectrum to fix the masses of the remaining particles, and then vary the ijn mass by hand. The 
solid lines from top to bottom in both (a) and (b) correspond to A em = 0, 0.3, 0.1. The dotted line 
is the nominal UED case from Fig. ^. 

line is the nominal UED case from Fig. ^, for which A^ m = 0.01. As expected, all lines 
follow the general trend of Fig. ||. In accord with the observations of Ref. [6] , we see that 
Iri coannihilations increase the prediction for £lh 2 . Such a behavior may seem peculiar at 
first sight, since in supersymmetry one finds the opposite phenomenon — coannihilations 
with sleptons tend to reduce the SUSY WIMP relic density. The difference between the 
two cases can be intuitively understood as follows. In SUSY, the cross-section for the 
main annihilation channel (xiXi ~^ ff) 1S helicity suppressed, but the coannihilation 
processes are not. Adding coannihilations therefore can only increase the effective cross- 
section ( |3.9| ) and correspondingly decrease Qh 2 . In contrast, in UED the main annihilation 
channel (7171 — > //) is already of normal strength. The effect of coannihilations can be 
easily guessed only if the additional processes have either much weaker or much stronger 
interactions. In the case of £ri, however, the additional processes are of the same order 
(both 71 and £bx have hypercharge interactions only) and the sign of the coannihilation 
effect depends on the detailed balance of numerical factors, which will be illustrated in 
Fig. H and discussed in more detail below. 

The spread in the lines in Fig. || is indicative of the importance of the coannihilations. 
Comparing Fig. ^|a and Fig. |||b, we see that in the case of three generations, the effects are 
magnified correspondingly. A similar conclusion was reached in Ref. [6]. 

Notice the peculiar ordering of the lines corresponding to different . With respect 
to variations of Ai R1 , the maximum possible value of Qh 2 is obtained for A^ m — ► 0, where 
the effect of coannihilations is maximal. Then, as we increase the mass splitting between l^x 
and 71, at first O/1 2 decreases (see the sequence of A^ m = 0, A^ m = 0.01 and A^ m = 0.1) 
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Figure 5: Plots of various quantities entering the LKP relic density computation, as a function 
of the mass splitting A eR1 between the LKP and the SU (2) \y -singlet KK electron, for = 500 
GeV in MUED. (a) Relic density, (b) a e ff(xp), (c) g e ff(xp) and (d) a e ff(xF)geff( x F)- 



but then starts increasing again and the Qh 2 values that we get for A^™ = 0.3 are slightly 
larger than those for Ai m =0.1. This behavior can be seen more clearly from Fig. ^a, 
where we vary the mass of the 5'C/(2)iy-singlet KK electron e_Ri and plot Qh 2 versus A ej?1 
for a fixed Rr 1 = 500 GeV. 

The interesting behavior of Qh 2 exhibited in Fig. [|a can be understood in terms of 
the mi R1 dependence of the effective annihilation cross-section ( |3.9| ) which is dominated 
by its a-term ( 3.17 ). Both a e fj and a e fj are functions of x, but for the purposes of our 
discussion here it is sufficient to concentrate on the fixed value x = xf which dominates 
the integrals ( |3,13j ) and ( |3.14 ). We plot a e ff(xp) as a function of A ejn in Fig. ||b. We see 
that a e ff(xF) exhibits exactly the opposite dependence to Qh 2 , and in particular, has an 
analogous local extremum at A eR1 ~ 0.1. Therefore, in order to understand qualitatively 
the behavior of flh 2 , we only need to concentrate on a e ff(xp). 

Let us start with the large A em region in Fig. ||b. The 5f7(2)w-singlet KK electron 
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Figure 6: The same as Fig. ^ but illustrating the effects of varying the SU (2) w -doublet KK electron 
mass. From top to bottom, the solid lines show Q,h 2 as a function of R~ l , for A eL1 = 0.01, 0.001, 0. 
The dotted line is the nominal UED case from Fig. ^. 



em is then too heavy to participate in any relevant coannihilation processes. The effective 
cross-section ( |3.9| ) then receives no contributions from processes with em- Nevertheless, 
the mass of em enters a e ff through the cross-section for the process 7171 — > e + e~ (see 
eqs. Q4.3D and Q4.5D ). Then as we lower m eR1 , 17(7171 — ► e + e~) is increased and this leads 
to a corresponding increase in a e fj as seen in Fig. |5|b. This trend continues down to 
A eR1 ~ 0.1, where coannihilations with em start becoming relevant. This can be seen in 
Fig. where we plot g e //(xp) as a function of A em . From its defining equation ( |3.10D we 
see that g e ff(x) starts to deviate from a constant only when the exponential terms (which 
signal the turning on of coannihilations) become non-negligible. The exponential terms are 
all positive and increase g e ff- At the same time, there are new cross-section terms entering 
the sum for o~ e ff, so we expect the numerator in ([3l]) to increase as well. This is confirmed 
in Fig. ||d, where we plot the numerator of (|3.9|) simply as a e ff(xF)g'jff{xF)- From Figs. ||c 
and ||d we see that both the numerator and the denominator of fl3.17|) increase at low A em , 
and so it is a priori unclear how their ratio will behave with A eR1 . In this particular case, 
g e ff wins, and a e ff(xp) is effectively decreased as a result of turning on the coannihilations 
with em- This feature was also observed in Ref. [6]. 

We are now in position to repeat the same analysis, but for the case of the SU(2)yy- 
doublet KK leptons £lx- In Fig. ||, in complete analogy to Fig. ||, we illustrate the effects on 
the relic density from varying the SU (2) ^-doublet KK electron mass. From top to bottom, 
the solid lines show f2/i 2 as a function of R , for A eL1 = 0.01,0.001,0. The dotted line is 
again the MUED reference model. We see that the case of SU(2) ^-doublet KK leptons is 
different. Unlike lm, they have weak interactions, and the extra terms which they bring 
into the sum ([T^) are larger than the main annihilation channel. The increase in g e jf is 
similar as before. As a result, this time the increase in the numerator of ( |3.S| ) wins, and the 
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Figure 7: The same as Figs. \^(b) and ^(b), but for the case of KK quarks. Each solid line is 
labeled by the value of (a) A 9JU or (b) A qL1 used. The dotted line is the nominal UED case from 
Fig. | 

net effect is to increase the effective annihilation cross-section. This leads to a reduction in 
the predicted value for the relic density, as evidenced from Fig. ||. Notice how the decrease 
in Qh 2 is monotonic with A eL1 . 

Another difference between £r% and £li coannihilations is revealed by comparing the 
case of 1 generation (panels (a) in Figs. || and |6|) and 3 generations (panels (b) in Figs. || 
and |). We see that for SU{2) W singlets, the coannihilations are more prominent for the 
case of 3 generations, while for SU(2)w doublets, it is the opposite. This is due to the 
different number of degrees of freedom contributed to g e ff in each case, which shifts the 
delicate balance between the numerator and denominator of (|3.9| ), as discussed above. 

5.2 Effects due to coannihilations with KK quarks and KK gluons 

We will now consider coannihilation effects with colored KK particles (KK quarks and 
KK gluons). Since they couple strongly, we expect on general grounds that the effective 
annihilation cross-sections will be enhanced, and the preferred range of the LKP mass 
will correspondingly be shifted higher. These expectations are confirmed by our explicit 
calculation whose results are shown in Figs. ^ and |8[ In Fig. [7] we show the effects on 
the relic density from varying the masses of all three generations of (a) jf>i7(2)w-singlet 
KK quarks and (b) «S?7(2)w-doublet KK quarks. The solid lines show Qh 2 as a function 
of R , and are labeled by the corresponding value of A qm or A qL1 used. As before, the 
dotted line is the MUED reference model. Comparing the results in Figs. 0a and ^b, we 
find that the coannihilations with qm and qn have very similar effects, as they are both 
dominated by the strong interactions, which are the same for qjn and qn- Fig. || shows 
the analogous result for the case of varying the KK gluon mass, where the labels now show 
the values of A ffl . There is a noticeable distortion of the lines around Z? -1 ~ 2300 GeV, 
which is due to the change in g* (see Fig. ^a). 
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Figure 8: The same as Fig. 0, but for the case of varying the KK gluon mass. The lines are labeled 
by the value of A 91 . The dotted line is the nominal UED case from Fig. || 



From Figs. [7| and |8] we see that in non- minimal UED models where the colored KK 
modes happen to exhibit some sort of degeneracy with the LKP, multi-TeV values for 
m 7l are in principle possible. From that point of view, unfortunately, there is no "no-lose" 
theorem for the LHC or ILC regarding a potential absolute upper bound on the LKP mass. 

5.3 Effects due to coannihilations with electroweak KK bosons 

We finally show our coannihilation results for the case of electroweak KK gauge bosons (W® 
and W^) and KK Higgs bosons (H®, G\ and Gf). The results are displayed in Figs. ||a 
and|9|b, correspondingly. Due to the SU(2)w symmetry, all three n = 1 KK IV-bosons are 
very degenerate, and we have assumed a common parameter for all three. Similarly, 
the masses of the n = 1 KK Higgs bosons differ only by electroweak symmetry breaking 
effects, which we neglect throughout the calculation. We have therefore assumed a common 
parameter Ajj for them as well. 

Since both the electroweak KK gauge bosons and the KK Higgs bosons have weak 
interactions, we expect the results to be similar to the case of St7(2)vi/-doublet leptons in 
the sense that coannihilations would lower the predictions for O/i 2 . This is confirmed by 
Fig. ||. We observe that the effects from the KK H^-bosons are actually quite significant, 
and can push the preferred LKP mass as high as 1.4 TeV. 



6. Summary and Conclusions 

In this paper we revisited the calculation of the LKP relic density in the scenario of Uni- 
versal Extra Dimensions. We extended the analysis of Ref. [6] to include all coannihilation 
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Figure 9: The same as Fig. [|, but illustrating the effect of varying simultaneously the masses of 
all (a) SU (2) KK gauge bosons and (b) KK Higgs bosons. In (a) the lines are labeled by the value 
of A Wl , while in (b) the values of A H are (from top to bottom) A H = 0.05,0.01,0.001,0. The 
dotted line is the nominal UED case from Fig. ||. 



processes involving n = 1 KK partners. This allowed us to predict reliably the preferred 
mass range for the KK dark matter particle in the Minimal UED model. We found that in 
order to account for all of the dark matter in the universe, the mass of 71 should be within 
500 — 600 GeV, which is somewhat lower than the range found in [6]. This is due to a 
combination of several factors. Among the effects which caused our prediction for Q,h 2 to 
go up are the following: we used a lower value of g* , we kept the individual KK masses in 
our formulas, and we accounted for the relativistic correction ( 3.1S| ). On the other hand, 
as we saw in Section [5], including the effect of coannihilations with KK particles other than 
/f>{7(2)vy-singlet KK leptons, always has the effect of lowering the predicted Qh 2 . Finally, 
the cosmologically preferred range for itself has shifted lower since the publication 
of [6]. 

The lower range of preferred values for R~ l is good news for collider and astroparticle 
searches for KK dark matter. It should be kept in mind that it is quite plausible, and in fact 
very likely, that the dark matter is made up of not one but several different components, 
in which case the LKP could be even lighter. We should mention that several collider 
studies [4,25,32,33] have already used an MUED benchmark point with = 500 GeV, 
a choice which we now see also happens to be relevant for cosmology. 

In Section || we also investigated how each class of n = 1 KK partners impacts the 
KK relic density. We summarize the observed trends in Fig. 10, where we fix fi/1 2 = 
0.1 and then show the required for any given Aj, for each class of KK particles. 
We show variations of the masses of one (red dotted) or three (red solid) generations of 
5C/(2)^y-singlet KK leptons; three generations of SU(2) ^-doublet leptons (magenta); three 
generations of 5C/(2)vi/-singlet quarks (blue) (the result for three generations of SU(2)\y- 
doublet quarks is almost identical); KK gluons (cyan) and electroweak KK gauge bosons 
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Figure 10: The change in the cosmologically preferred value for as a result of varying the 
different KK masses away from their nominal MUED values. Along each line, the LKP relic density 
is fl x h 2 = 0.1. To draw the lines, we hrst fix the MUED spectrum, and then vary the corresponding 
KK mass and plot the value of which is required to give tt x h 2 = 0.1. We show variations of 
the masses of one (red dotted) or three (red solid) generations of SU (2) w singlet KK leptons; three 
generations of SU (2) \y -doublet leptons (magenta); three generations of SU(2)w singlet quarks 
(blue) (the result for three generations of SU (2) \y -doublet quarks is almost identical); KK gluons 
(cyan) and electroweak KK gauge bosons (green). The circle on each line denotes the MUED values 
of A and RT 1 . 

(green). The circle on each line denotes the MUED values of A and R^ 1 . 

Fig. |l(] summarizes our results from Section ||. It also provides a quick reference 
guide for the expected variations in the predicted value of Qh 2 as we move away from the 
Minimal UED model. For example, it is clear that unlike the case of coannihilations with 
£rx, which was considered in [6], coannihilations with all other KK particles will lower the 
prediction for Qh 2 and correspondingly increase the preferred range of R . This is due to 
the larger couplings of those particles. Fig. |l^ can also be used to quantitatively estimate 
the variations in the preferred value of R -1 in non-minimal models. 

On a final note, in the non-minimal UED model, other neutral KK particles such as Z\ 
can also be dark matter candidates. On dimensional grounds, the relic density is inversely 
proportional to the square of the LKP mass, 



Due to the larger coupling #2 of the SU (2)w gauge interactions, we expect the upper bound 



nh 2 



Wi 2 




(6.1) 



(6.2) 
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on mz 1 - l consistent with WMAP, to be larger than the bound on m 7l roughly by a factor 
of g^/gi ~ 3. However, in the Z\ LKP case, SU(2)w symmetry implies that the charged 
W\ states are almost degenerate with Z\, and therefore coannihilations with will be 
very important and will need to be considered. The analysis of the cases of Z\ and H\ 
LKP and their detection prospects is currently in progress [71]. The results presented in 
this paper are also relevant for the case of KK graviton superwimps [72-74], whose relic 
density is still determined by the freeze-out of the next-to-lightest KK particle. 

In conclusion, dark matter candidates from theories with extra dimensions should be 
considered on an equal footing with more conventional candidates such as SUSY dark mat- 
ter or axions. The framework of Universal Extra Dimensions provides a useful playground 
for gaining some experience about the signals one could expect from extra dimensional dark 
matter. If extra dimensions have indeed something to do with the dark matter problem, 
the explicit realization of that idea may look quite differently (see for example [75-77]), 
especially if one wants to resolve the radion stabilization problem [78,79]. Nevertheless, we 
believe that the methods and insight we developed in this paper will prove useful in more 
general contexts. 

Note added: During the completion of this work we became aware of an analogous 
calculation done independently by another group [12]. We have compared extensively the 
formulas for the annihilation cross-sections involving KK quarks, KK leptons and KK gauge 
bosons, in the limit of degenerate KK masses, as listed in the Appendix. In all considered 
cases we found perfect agreement. 
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A. Annihilation Cross-Sections 

In this section, we summarize the annihilation cross sections of any pair of n = 1 KK 
particles into SM fields in the limit of no electroweak symmetry breaking, as in [6]. In 
order to render the formulas manageable for publication, in this Appendix we list our 
results in the limit of equal KK masses. However, in our numerical calculation, we kept 
different masses for all KK particles, which often leads to enormously complicated analytical 
expressions. We also assume all SM particles to be massless, since we are working in the 
limit where we neglect electroweak symmetry breaking (EWSB) effects of order vR, where 
v is the Higgs vacuum expectation value in the SM. All cross-sections are calculated at 
tree level. All vertices satisfy KK-number conservation and KK-parity since KK-number 
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violating interactions are only induced at the loop level [5] . Some of the cross-sections have 
already appeared in [6] and we find perfect agreement with those results. We define a few 
constants below which are commonly used in our formulas for the cross-sections. 

91 = — , (A.l) 

92 = -, (A.2) 
gz = ^— , (A.3) 



/3 = yi-^ J (A.4) 

L = lo s (ir!) = " 2tanh " 1/3 • (A - 5) 



Here g±, g^ and are the gauge couplings of U(l)y, SU(2)w and SU(3). (3 is the velocity 
of the incoming KK particle in the annihilation process. Notice that L is negative since 
< /3 < 1. m is the KK mass which for the purposes of this appendix is the same for all 
KK particles, e is the electric charge and c w and s w are the cosine and sine of the Weinberg 
angle in the SM. Table [l| provides a quick reference guide for the different process types. 
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Table 1: A guide to the formulas in the Appendix. Each box in the table corresponds to a 
particular type of an initial state. The entry points to the section in the Appendix where the 
corresponding annihilation cross-sections can be found. Here "gauge bosons" include EW KK 
gauge bosons (W* , Z\ and 71) and the KK gluon (gi). "Higgses" stands for the KK Higgs (Hi) 
and KK Goldstone bosons (Gf and Gi). Leptons contain both SU (2) w singlet KK leptons (£ri) 
and ££/(2)w/-doublet KK leptons (£1,1 and v^). Quarks include both SU (2) ^-doublet KK quarks 
(qui) an d SU (2) w -singlet KK quarks (<Zli). 



A.l Leptons 

Coannihilations with S'L r (2)vi/-singlet KK leptons £r\ are important since they are expected 
to be the next-to-lightest KK particles in the Minimal UED model [5,6]. For fermion final 
states with / 7^ £, the cross-section is 

N c gfY f 2 (Yj +Y?)(s + 2m 2 ) 
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where N c is 3 for quarks and 1 for leptons. For cases with the same lepton flavor in the 
initial and final state, we have 



+ 



+ 



gjY e A Ri (5(3s + 2(2s + 3m 2 )L) 

32vr/? 2 s 2 

9l Y e 4 m (/?( 4s + 9 ™ 2 ) + 8m2L ) 
647rm 2 /3 2 s 



(A.7) 



t ft , D±B 



t e ±\ 



24vr/3s 2 

gfY e 4 (-m 2 (4s - 5m 2 )L - f3s{2s - m 2 )) 
32vr/3 2 s 2 m 2 

... gtY^s-Zm 2 ) 



(A.8) 
(A.9) 
(A.10) 



647r/3 2 sm 2 

where I and are the leptons from different families. For the remaining final states we get 

^Y/Y 2 (s + 2m 2 ) 



a(£+£ 



gfY e 4 (2(s 2 + 4m 2 s - 8m 4 ) - /3s(s + 4m 2 )) 



(A.ll) 
(A.12) 



Our results, ( |A.6| - [A.12"|) , exactly agree with (C.l) - (C.8) from [6]. 

The cross-sections among left handed fermions are somewhat complicated since they 
involve SU(2)]y gauge bosons as well as U{1)y gauge bosons. For KK neutrinos we find 



cr(mm -> ff) 



N c g 2 z (gj + g 2 R )(s + 2m 2 
2Air(3s 2 

g 2 g 2 (s + 2m 2 ) 



1 



o{vnVti -> ZZ) = 
o{yi\vi\ — > W + W 



48tt/3s 2 

#f ((8m 4 - s 2 - Am 2 s)L - (3s{s + 4m 2 )) 
8vr/3 2 s 3 

_ 5g|( g + 2m 2 ) gl((3s-2m 2 L) 

' 96vr/3s 2 32vr/? 2 s 2 

g%(f3(s + 4m 2 ) + (s + 2m 2 )L) 
327T/3V ' 

g%{(3s{2s - m 2 ) - m 2 (4s - 5m 2 )L) 



(A.13) 

(A.14) 
(A.15) 



327r/3 2 s 2 m 2 
5 |(4s - 3m 2 ) 



647r/3sm 



2 ' 



(A.16) 

(A.17) 
(A.18) 
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, - g|(/?(4 g + 9m 2 ) + 8m 2 L) 

<y\yi\vi>\ -> tw) = — ^— ^ — 2 ' (A. 19) 



„(„ n ^e~P>+\ 9tms + 9m 2 )+8m 2 L) 

o\yi\vv\ -+t t ) = 2 — 5 > (A.20) 

a(^a - * o = 

g 2 z (gl+f R )(s + 2m 2 ) gj(/?(4 g + 9m 2 )+8m 2 L) 

24vr/?s 2 647rm 2 ^ 2 s ' 1 ' ' 

Here = j^- (T 3 — Q js 2 ) , §l = 9z for neutrinos and gz, = g2/V% for charged leptons. 
94> = T^c~ ~~ Q<t> s w) w ^ n Q<p = 1 fo r the upper entry in the Higgs doublet and Q<p = 
for the lower entry. Since we ignore EWSB, all gauge bosons have transverse polarizations 
only. 4> represents either a charged Higgs boson (4> u , isospin i) or a neutral Higgs boson 
(<j> d , isospin 

The previous results allow us to immediately obtain 



v{veim -»■ <P4>*) 






<y{yvJ>i\ -»■ ZZ) 


= ^(4i^Zi ■ 


-> ZZ + Z7 + 77) 


a( m u a W+WH 


= ^(4i^Zi ■ 


-»• W + W~) , 




= ^iAi 


, 




= o-(^a4i - 


- £ ± / ± ) , 



(A.22) 



For at least one charged KK lepton in the initial state we get 



^ZAl ^ 0r = 327r/? 2 s 2 

ffl(/3(4 g + 9m 2 ) + 8m 2 L) g 4 ( S + 2m 2 ) 

64vrm 2 /3 2 s 24vr/3s 2 ' 1 ' ' 

a(£ L1 u L1 -.//)= 967r/3s2 - (A-25) 

- «W = , (A-26) 

ff r/",7 -.IV-Rl gfg 2 2 ((8m 4 - S 2 - 4m 2 s)L - + 4m 2 )) 
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a(£ L1 u L1 W 3 ) = 48 ^ 2/3 + 327Tg2(32 



gj{(3(s + 4m 2 ) + (s + 2m 2 )L) g 4 m 2 L 

64vr S 2 /3 + i6vr S 2 ' 1 8j 



g lg2 3 (5/? g + 2(2 g + 3m 2 )L) 
a(£ L1 u L1 - / i*) = 64 ^ 2s2(24 _ 1} 



+ 



g 2 g 2 (/3(4 S + 9m 2 ) + 8m 2 L) g|(g + 2m 2 ) 

64vrm 2 /? 2 s(24 - l) 2 + 96tt/?s 2 ' 1 ' 



gig^(-2m 2 (4:S - 5m 2 )L + m 2 Ps) 



64vr/? 2 s 2 m 2 (24 - 1) 
M4s-3m 2 ) f g 2 *? 2 , + ( A 



64vr/3 2 s 2 m 2 V( 2 4-!) 2 4 

^ } = 256^/3 ' (A " 31) 
^„ f -w,,^ g 2 ffl(4s-3m 2 ) 

^i/xi - ^ ) = 64Ws/3(2s 2 _ 1)2 > (A.32) 

p , ff 2 g 2 (/3(4s + 9m 2 )+8m 2 L) 

" W£1 ^ ^ } = 64 7 rm 2 S /? 2 (24 - I) 2 ' (A " 33) 

a ^ _x gt(/3(9m 2 + 4 S ) + 8m 2 L) 

- w) = 256 vrm 2 s/ 3 2 ' (A " 34) 



where g 2 = g 2 YyY£ fl + g^TjT^. The above cross-sections, ( A.13 - A. 34 ) are consistent with 
(B.48) - (B.62) and (B.71) - (B.74) in [6]. 

For one S't r (2)vi/-singlet KK lepton and one S'C r (2)vi/-doublet KK lepton we get 

4y2 y2 

a{£ R1 i L1 ^££) = (8m 2 L + /?(9m 2 + 4s)) , (A.35) 

oVrlIli - U) = ( 4s - 3m ') • (A. 36) 
These two formulas have the same structure as <^(£%\£ ,: ri — * £ ± £' ± ) and^/m l ± ^' T ). 
A. 2 Gauge bosons 

The self-annihilation cross-sections of 71 are 

q4(Y 4 + y 4 ) 

*(TL7i-//) = 72^2^2 7 (-5s(2m 2 + s)L-7s/?) , (A.37) 

^(7171 - = ^ • (A.38) 
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These two cross-sections are identical to (A.44) and (A. 47) in [6]. For Z\ self-annihilation 
into fermions and Higgs bosons, 

cr^Z, - //) = ("5(2m 2 + s)L - 7s0) , (A.39) 

- tStv (a ' 40) 

The cross-section for the above two processes are obtained from (7(7171 — »■ //) and 
(7(7171 — ► 00*) by replacing g{Y with 52/2, which corresponds to the Z couplings to 
SM fermions and Higgs bosons. For the coannihilations of SU(2)w KK bosons into SM 
gauge bosons, we get 



aiZiZx ^ W + W~) = 9 l , n0 (l2m 2 (s - 2m 2 )L + s/3(12m 4 + 3sm 2 + 4s 2 )) 

18irm z s d p z 



(A.41) 



a(W+W+ -► VF+VF+) = (I2m 4 (s - 2m 2 )L + s/3(12m 4 + 3sm 2 + 4s 2 )) 

36irm 2 p z 



e 4 

^W^Wf 77) = 3gWg3 ^ 2 (I2m 4 (s - 2m 2 )L + s/?(12m 4 + 3sm 2 + 4s 2 )) , (A.43) 

222 

a(W+Wr 7 Z) = 92G C ™ n0 (I2m 4 (s - 2m 2 )L + s/3(12m 4 + 3sm 2 + 4s 2 )) , (A.44) 
I8irm z s 6 p z ' 

4 4 

afWi+Wr ZZ) = 92 ^°, n0 (I2m 4 (s - 2m 2 )L + s(3{12m 4 + 3sm 2 + 4s 2 )) . (A.45) 

36-^77^ s 6 f3 z 

We see that the above five cross-sections contain similar expressions up to overall factors 
due to the gauge structure of SU (2)w x U(l)y ■ For annihilation into other final states, 
we have 

*(W+Wr^ff) = ((12m 2 + 5s)L + 2/3(4m 2 + 5s)) , (A.46) 

a 4 

a(W+Wi^W + W-) = 187r ^2 g2/3 2 ( 2 ^ 2 (3m 2 + 2s)L + /3(llm 4 + 5sm 2 + 2s 2 )) , 

(A.47) 

«-« = ^. (A.48) 

These three cross-sections are different since they involve s-channel Z diagrams. For 71 Zi 
and iiWi into fermions we can recycle (7(7171 — > //) and obtain 

-N c g 2 g 2 Y 2 

a(7i Zi - //) = 288 ^/ (5(2m 2 + s)L + 7s/?) , (A.49) 
^(TlWT ^ //') = 1447rg2 y (5(2m 2 + s)L + 75/3) . (A.50) 



(A.42) 
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For the annihilation of two different KK gauge bosons into Higgs bosons we have 

^ z ^^ = Mp- (A - 5i > 

*(nwr - Mi) = J$ ■ (A - 52) 

aiZtWr - <fa« = , (A.53) 

which can be obtained from (7(7171 — ► The cross-section for Z\W^ into fermions 

^(Zi Wf //') = 5 ~ 67r ^ 2 ((14m 2 + 5s)L + /3(16m 2 + 13s)) , (A.54) 

has a different structure compared to other fermion final states due to an s-channel W 
diagram. The cross-sections for Z\W{ into gauge boson final states 

4 2 

a(ZiW{ — > ZW~) = 187r ^2g 2/3 2 ( 2 ™ 2 (3m 2 + 2s)L + /3(llm 4 + 5sm 2 + 2s 2 )) , (A.55) 

2 2 

<t(Zi Wf = - 6 f% _„ (2m 2 (3m 2 + 2s)L + p(llm 4 + 5sm 2 + 2s 2 )) , (A.56) 

187rm z s z p z 

can be obtained from aiW^W^ — ► H^ + M^ _ ). For KK gluons we get 

o 4 

vigigi^gg) = . I (8m 2 (s 2 + 3sm 2 - 3m 4 )L + sf5(Um 2 + 13sm 2 + 8s 2 )) , 

(A.57) 

<K<7i<7i -<?«) = ^^2 (2(20s + 49m 2 )L + /3(72m 2 + 83s)) , (A.58) 

for which there are no analogous processes. The cross-sections associated with one gluon 
and one electroweak gauge bosons in the initial state 

^171^99) = ^^^(-5(2m 2 + s)L-7s/3) , (A.59) 
a(SiZi-<ra) = (-5(2m 2 + s)L-7s/3) , (A.60) 

2 2 

a(g 1 W 1 - -+ qq>) = 2 ™* 202 ("5(2m 2 + s)L - 7s/?) , (AM) 

are obtained from (7(7171 — > //), (7(71^1 — > //) and (7(71 W-f — > //') by simple coupling 
replacements and accounting for the additional color factors. 
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A. 3 Fermions and Gauge Bosons 

Note that in UED, the Dirac KK fermions are constructed out of two Weyl fermions with 
the same SU(2)w x U (l)y quantum numbers while a Dirac fermion in the Standard Model 
is made up of two Weyl fermions of different SU(2)w x U(l)y quantum numbers. Therefore 
the couplings of KK fermions to zero mode gauge bosons are vector-like. This difference 
shows up in processes involving gauge-boson couplings with fermions. For the vertices 
which involve n = 1 gauge bosons, we need one KK fermion and one SM fermion in order 
to conserve KK number. In this case, there is always a projection operator associated with 
the subscript (L/R) of the KK fermion. The annihilation cross-sections with SU(2)\y- 
singlet KK fermions and SU(2)w KK gauge bosons {Z\ and W^) are zero: 

aifmZi - SM) = , 
otfiaWf ^SM)=0. (A.62) 

The cross-section for 577(2) ^-singlet leptons and 71 is 



For a KK quark and 71 we have 

viQili^qg) = ^Jf (-HL- )) . (A.G1) 



This cross-section is basically the same as 17(71 — ► BqI^), up to a group factor. In 
f(<Zi7i qg), the vertex associated with the SM gluon g contains a Gell-Mann matrix t®-. 
In the squared matrix element, we then get [80] 

E^E^ = ^4 x3 = ^ (A. 65) 

Similarly, for the 577(2) v^-doublet quarks with Zi, we get 

o{q L1 Z x gq) = ^JgL ( _ 6L _ p) , (A.66) 

by replacing the g\ coupling in a(qi^i — > qg) with 52/2. For W^, one should use the 
coupling g2/V2 instead: 

o{q^W, - W) = (-6L - f>) • (A.67) 

For the 5C/(2)jy-doublet KK leptons and 71 or Z\, we get 
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- "<W) = (-6L - ffl . (A.70) 

The last 7 cross-sections have a similar structure since they all have s- and ^-channel 

diagrams only. For the cross-sections associated with 5C/(2)^-doublet KK leptons and 
electroweak KK gauge bosons into other final states, we have 

o{l LX Z x ^v t W) = ^-^2^2 (26m 2 L + /3(23m 2 + 32s)) , (A.71) 

a(e LlWl -+ Ml + Z)) = 768 J s p cl K(324 - 6)L 

+/3m(24c 2 - 1) + 32s/3c 2 ) , (A.72) 

a(fLiWf - IlW~) = l927T 9 l 2s(32 (-3m 2 L + 4*0) , (A.73) 

^li^+^W 1 ") = 3847r S^ (l6m 2 L + /?(llm 2 + 8 g )) , (A.74) 

ff^LiWT -► W) = -» »iW + ) , (A.75) 

a(^Li^i + £lW + ) = a{v tl W{ -► i/^") . (A.76) 

For KK gluon - KK quark annihilation, we obtain 

cr(gm^gq) = (24m 2 L + /3(25m 2 + 36s)) . (A.77) 

A. 4 Quarks 

The annihilation cross-section of two KK quarks into SM quarks of different flavor 

( - '-'\ gj{s + 2m 2 ) 

*(Qm - * 9 ) = Uns2p (A.78) 



can be obtained from c(^ri^ri ~* ff) by multiplying with the following group factor [80] 
\\tr (t a t b ) tr (t b t a ) = l -C{r)H ab 5 ab = ± ■ \ • 8 = I. (A.79) 



Here C(r) = 1/2 is the quadratic Casimir operator for the fundamental representation of 
SU(3). KK quark annihilation into same flavor SM quarks is given by 

a(q iqi ^qq) = 432^2^2 {^ 2 {As - 5m 2 )L + (6s - 5m 2 )s(3) . (A.80) 

In this process, there are three terms in the squared matrix element since we have both 
t- and u-channel diagrams. This process also has an analogy with c^m^m — ► £ ± £ ± ). 
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However, each term gets a different group factor. The squares of the t- and u-channel 
diagrams get the same factor of 2/9 but for the cross term we obtain [80] 

l ^j\r (ViW) = l - (c 2 {2) - \C 2 {G)^ tr(t a t a ) 



1/4 3\ 4 _ 2 _ 2 
9 \3 ~ 2 ) 3 ~ ~9 ~ ~27 



(A.81) 



For qiqi annihilation into same flavor SM quarks we get 



o 4 

<j{qiqi^qq) = ^^^2 (4m 2 (4s - 3m 2 )L + /3(32m 4 + 33sm 2 + 12s 2 )) , (A.82) 

which can also be obtained using the analogy to o^f^f^ — > £ + £~) and taking into account 
group factors. For the final state with gluons we have 

<r(qiqi^99) = (4(m 4 + 4sm 2 + s 2 )L + s/3(31m 2 + 7s)) . (A.83) 

For different quark flavors in the initial state, we have 

i i i\ <7t(4s — 3m 2 ) 

^ ^ = W ' (A - 84) 

- tfO = 288 ^ 2g/32 (8m 2 L + /3(9m 2 + 4s)) . (A.85) 

The above two cross-sections can be obtained from cr(£^ 1 ^'^ 1 — > ^ =t £ /= ' : ) and cr(£^ 1 £'^ 1 — > 
^ =t £ /=F ), correspondingly, by multiplying with the group factor 2/9. The remaining cross- 
sections are 

= o(uRid L i — »■ ueZ) , (A. 86) 

= a(u R1 u L1 -»■ u«) , 

and 

0"(<?i<?i — >■ w') = cr(uRid,Ri -»■ ltd) , 

= a(u R id L1 -> W) , (A.87) 
= a(u R1 u L1 -> uu) . 

A. 5 Quarks and Leptons 

The cross-sections listed below are mediated by t- or u-channel diagrams with KK gauge 
bosons. For one KK lepton and one KK quark in the initial state, we get 

- fa) = ^g^g i , (A.88) 
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a(v lUL1 - vu) - 102Ws/3 > (A. 89) 

m a - 3m 2 ) 

- ,d) = • (A.90) 

These three cross-sections can be obtained from cr(£^ 1 £'% 1 — > £ ± £ /± ). For one KK lepton 
and one KK anti-quark in the initial state, we have 

a(v lUL1 - M) = 256 ^ 2s/32 (8m 2 L + /3(9m 2 + 4,)) , (A.91) 

a(£ L1 u L1 - te) = (4 ^ ^ 2 ; /3 ? 2)2 (8- 2 L + /3(9m 2 + 4,)) , (A.92) 

afrifci - vu) = { ' 9 ^l 2 + jF (8- 2 L + /?(9m 2 + 4,)) , (A.93) 

which can be obtained from ^{i^^W — ► ^ ± ^ /T ). If one of the particles in the initial state 
is an S'J7(2)^y-singlet fermion, only 71 can mediate the process and the cross-sections can 
be obtained from our previous results: 

ct(£riUli -> tu) = cr(£ L1 u R1 -> £u) = a(£ R1 u R i -> £tZ) , 

= <x(W hi -► = a(£ R1 £ L1 -► , (A.94) 



o(£r\ul\ — »■ ^u) = o(£l\ur\ — »■ £u) = o(£ R \u R \ — >■ £u) , 

= ct(^i4i *0 = ^(^meLi . (A.95) 

A. 6 Higgs Bosons 

The mass terms for the KK S'C/(2)vi/-singlets appear with the wrong sign in the fermion 
Lagrangian. For example, the mass term for the KK quarks leads to the following structure 
for the mass matrix at tree level 



(u Ln (x),u Rn (x)) 

R 171 \ I u Ln{x) 

\ m ~Tl) \ u Rn(x) 
The corresponding mass eigenstates u' Ln and u' Rn have mass 



(A.96) 



M - = V© 2+m2 - (A - 97) 

These mass eigenstates receive different radiative corrections which lift the degeneracy [5]. 
The interaction eigenstates are related to the mass eigenstates by 




cos a 75 sin a 
sin a —75 cos a 




(A.98) 
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where a is the mixing angle between SU (2)w-singlet and SU (2) ^-doublet fermions defined 
by 

Tfi 

tan2a = — -. (A.99) 
n/R 

This mixing is very small except for the top quark. However, even with a ~ the effect 
of the rotation ( A.98p is present in the Yukawa couplings through the redefinition UR n — > 
— 7 5 UR n . It does not affect the gauge-fermion couplings. We use the following notation for 
KK Higgs bosons, 

y h?^ Gi ) ■ (a.ioo) 

We keep only the top- Yukawa coupling and we also keep the Higgs self-coupling assuming 
rrih = 120 GeV. Below we list the cross-sections associated with two KK Higgs bosons in 
the initial state. 

a(H x H x ^ G+G~) = ^— j— (8e 2 m 2 X 2 s 2 w L + (3{{2s + m 2 )e 4 

+4\m 2 s 2 w e 2 + 4A 2 m 2 4» , (A.101) 
9A 2 

a(H x H x -HH) = —, (A.102) 

a(H lHl _> GG) = l2Sn J s p sici (8e 2 m 2 Xs 2 w c 2 w L + 0{(2s + m 2 )e 4 

+4Am 2 e 2 s 2 + 4A 2 m 2 4ct}) , (A.103) 
a(H x H x ^ZZ) = M7T ^ p2c4 {sP(s + 4m 2 ) + 4m 2 (s - 2m 2 )L) , (A.104) 

a(H x H x -> W + W~) = on g A 9 a + 4m 2 ) + 4m 2 (s - 2m 2 )L) , (A.105) 

3w 4 

a(H x H x ^tt) = - U ' 2 (-(s + 2m 2 )L-2s[3) , (A.106) 

<t(G+G+ - G+G+) = — * 44 (-16e 2 m 2 A S 2 4L 

1287rm z sp z s* c* v 

+ /5{(2s + m 2 )e 4 -8Am 2 e 2 s 2 c 2 y + 16A 2 m 2 s 4 y ct}) , (A.107) 



a(G+G^ - «) = 11527rg2/?2g4c 4 {72sic 2 w y 2 (Ssc 2 w y 2 - Am 2 e 2 )L - 432 S /?44y t 4 

+ s/3 3 (204 - 124 + 9)e 4 - UAsf3sic 2 w y 2 e 2 ) , (A. 108) 



a(GfG~[ -> 66) = 11527rg 2 / j2 s 4 c 4 (7244y|(-3 S 44y t 2 + e 2 m 2 (44 - 3))L 



-432s/344y 4 + s/3 3 (204 - 244 + 9)e 4 (A. 109) 

-36 S /342/ 2 e 2 (44-74 + 3)) , 
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a(G+G7 _> G+G-) = 1927rm2g2/ j 2g 4 c 4 (6m 2 e 2 S (e 2 + 4As 2 u c 2 )L - A8X 2 m 2 SS 4 w ci 



+ /3{e 4 (m 4 - 7sm 2 - 3s 2 ) - \2\rn 2 ss 2 w c 2 w e 2 }) , (A.110) 



(7(G + Gr _ = 1287rm a^ 4 (8 S > 2 e 2 AL + /?{(2 S + m 2 )e 4 

+4A4m 2 e 2 + 4A 2 m 2 s 4 y }) , (A.lll) 



(7(G + Gr ^ Gi7) = 4g7r J s ^ 2 (24m 2 c 2 yS L - /?{4(1 - 24) 2 m 4 

+ s(92s* - 1404 + 47) - 24s 2 c 4 u }) , (A.112) 

-A ( 1 _ Oo2 \4 



o{G\G\^ZZ) = ^—^L(4m 2 (s-2m 2 )L + s(3(s + 4m 2 )) , (A.113) 
<r(G+G7^ 7 ^) = 2*8*0* ^ 4m2(s ~ 2m2)L + + 4m2 ^ ' (A ' U4) 



e 4 



ct(G+G7^77) = ^3^2 (4m 2 (s - 2m 2 )L + + 4m 2 )) , (A.115) 
a(G+Gi ^W + W~) = 2i ^ 2(32 (6m 2 L + (3{s + 11m 2 )) , (A.116) 

<x(G+G7 - //) = ^ + &) , (A.117) 
where / represents leptons and first two generations of quarks and 

= -e 2 Q/ " 2^1(1 - 2s 2 )(Tf - Q f s 2 J. 
A number of cross-sections can be simply related: 



<t(GiGi -► G+G~) 


= (T^ifTi - 


-> G+G~) 


f(GiGi - 


► GG) 


= (T^ifTi - 


+ HH) , 


<t(GiGi — ► 




= <7(fflfTi - 


-GG) , 


<7(GiGi - 


> ZZ) 


= <7(fTifTi - 


+ ZZ) , 


(GiGi -► W 


h WH 






a(GiGi 


-> if) 


= - 




(j(G+G7 - 


► GG) 


= a(G+G7 ■ 


^HH) , 


<r(G+G7 - 


ffif) 


= a(G+G7 ■ 


^GG) , 



^a(HiHi - G+G- 



(A.118) 
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The rest of the cross-sections are 

,(H lGl EG) = 19W ' 2ffl44 (*»V.(.» - IXslcDL 

,4 v„™2 o„2\ 4 



-/3{(m 4 - 7sm 2 - 3s 2 )e 4 (A.119) 
+6Am 2 ss 2 c^e 4 — 12A 2 m 2 ss 4 c 4 }) , 



cr(HiGi - G+G-) = (24m 2 c 2 S L + /3{-4(l - 2 S 2 ) 2 m 4 

+W(-92st + 140s 2 -47) + 24s 2 ct}) , (A.120) 

aiHid ^ W + W~) = 96 J 2 3 ^ 2 (12m 2 (2m 2 + s)L + s/3(32m 2 + a)) , (A.121) 

v{H\G 1 - tt) = 2887rs2 ^ 2c2 g2 (-54y 2 (e 2 m 2 + c^4(s - 2m 2 )y 2 )L 

+ /3{-54 S y 4 C 2 y 4-27e 2 y 2 S + e 2 <7| S /3 2 (324-244 + 9)}) (A.122) 

a^G^ff) = ^{9l + 9%) , (A.123) 

a(^G+ - tfG+) = l927T J s2p2si (6m 2 e 2 s(e 2 - 2A S 2 )L - 12A 2 m 2 S 4 

- /3{(m 4 -7sm 2 -3s 2 )e 4 + 6Am 2 s4e 4 }) , (A.124) 



4 

487rm 2 s 2 /? 2 

,4 i /|f7„™2 0/t„2\4 



(7(^0+ - = - 9z 22Q2 (I2m 2 (l - 24)(2 - 3 2 )L 



-/3{(4m 4 + 47sm^-24s^)c 4 (A.125) 
-24(m 2 - s)ss 2 w c 2 w - 3s(m 2 + 4s)s 4 }) , 



aiH.Gt - ZW + ) = (6m 2 L{(l - 2s 2 )(4m 2 - 2ss 2 + s) + s} 

+ s/3{44(s + 11m 2 ) + (1 - 2s 2 w )(s + 32m 2 )}) , (A.126) 

2 2 

a(F!G+^ 7 ^ + ) = ^Jf^ (6m 2 L + /?(llm 2 + S )) , (A.127) 



<7(#iG+ - *6) = 64?rs2/32g4 (-64(2e 2 m 2 + SS 2 y 2 )y 2 L 

+ /3{ s/ 3 2 e 4 - 6ss 4 w y 2 e 2 - 12ss 4 w yt}) , (A. 128) 

a(H lG f - //') = ^gg , (A.129) 
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+ GG + ) 


= (7(^1 G+ - 


-> fTG" 


o{GiG+ - 




= a(H x G+ - 


-» GG 4 


a(G x Gf - 


ZV^+) 


= a( J ff 1 G+ - 


-» 


o{GiG+ - 


>7^+) 


= a( J ff 1 G+ - 




a(G x G+ 


- //') 


= a^Gf - 


- //"') 



(A.130) 



A. 7 Higgs Bosons and Gauge Bosons 

The cross-sections involving one KK Higgs boson and one KK gauge boson are given below. 
They are rather simple compared to the cross-sections from Section A.£. 



a{H l9l tt) 



ght 



4&nm 2 s 2 (3 2 
a{H 1 Z 1 -► ZH) 

a[H x Z x -» W~G + ) 



(2m 2 (s - m 2 )L + (2s - 5m 2 )s/3) 



4 



(L + 4/?) 



967rm 2 s 2 /3 2 



967T4S/3 2 

4m 2 L + s/3(4s + m 2 )) 



a{H X Z X ^iE) = 64 ^2 sf3 2 ( 2m L + ( 4fi - llm )P) 



2 \2 



a(GfZ x -» ZG+) 
^(G+Z! - 7 G 4 



9Q-kc 2 w s[3 2 



'w ' 
,2„2 



,2. ,2 



a^-tS) = (2m 2 L + (4, - 1 



<r(ffi7i -» tt) 



a(tf l7l - Zfl") 
o-(Gi7i -> W~G 

576irm 2 c 2 v s 2 f3 2 
a{G+ 11 ^ZG + ) = 

a(G+ 7 i - 7G + ) 
a(G|VF 1 + -» G + W + ) 



2 2 

#1#2 



96vrc 2 s/3 2 



(L + 4/3) , 



(L + 4/3) 



y 96vrs/? 2 

2m 2 (7s + 8m 2 )L + (4s - 43m 2 )s/3) 



96ircls/3 2 



e 2 5? 



(L + 4/3) , 



52 1 



96Trm 2 sf3 2 



2Ans(3 2 
(l2m 2 L + /3(6m 2 + 5s)) , 
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a(Hl wt G+Z) = (m»(2 - 4 - 24,)L 

-H{m 2 Hsi - 4 + 1) + (34 - 74 + 4)}) , (A.144) 

2 2 

<x(/W+ - G+ 7 ) = Q6 ^ 2gp2 {~2m 2 L + /3(4m 2 + 3s)) , (A.145) 

a(H 1 W+^HW + ) = -^(L + 4/?) , (A.146) 

o(G- lW ?^W + G-) = ^-|L(-2m 2 L + /3(4m 2 + 3 S )) , (A.147) 



a(G T Wt ZG) = 96 jl^ Kd24 - 154 + 4)L 
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+/?{(64-34 + l)m 2 + S (54-94 + 4)}) , (A.148) 

2 2 

a{G - W + _ lG) = QQ ^ 2gp2 (l2m 2 L + /3(6m 2 + 5s)) , (A.149) 

a (GiW+^tt) = 647r f 2g/3 2 (-^VtL + /?{e 2 m 2 + 2(4s - llm 2 )^ 2 }) , (A.150) 

where yt is the top quark Yukawa coupling. The cross-sections listed below are obtained 
from our previous calculations. For one KK Higgs boson and one KK gluon we get 

a(H l9l tt) = a{G l9l tt) , 

= a(G+ 9l - 46) . (A.151) 

For one KK Higgs boson and one 2q we have 

a(H 1 Z 1 Zif) = (t(GiZi Zfl") , 
<7(fTiZi -»• W^G+) = <t(GiZi -► W~G+) , 
= <r(Gf Zi -»• VF+G) , 
= <r(Gf Zi -► VF+#) , 

= (7(^+^1 -► W+G) , (A. 152) 

= (7(^+^1 -► W + #) , 
<7(fTiZi -► tt) = cr(GiZi -► tt) , 
(y{G\Z x -► i&) = <7(iWi + -► *6) , 

= <r(GiW l + -► *6) . 
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-> G+Z) 


= ct(GiW7 


-»• G+z; 


a{H{W^ - 




= <r(GiW l + 


-G+7) 


a{H{W^ -» 


FVF+) 


= a{G{W^ 


-► GVF+ 




-»• ZG) 


= <j(GiW? 


-► Zff) 


a(G^W+ 


-7G) 


= a{G^W+ 


^7#) 



For one KK Higgs boson and one 71 we obtain 

a(H l7l - Zff) = a(Gi 7 i - ZG) , 
<r(Gi7i -► VF"G + ) = a(iJi7i -»• , 
= a(G+ 7l - W+G) , 

= a(G+7i -► VF+#) , (A.153) 
<7(fTi7i tt) = a(Gi7i tt) , 
= a(G+ 7l - *6) . 

For one KK Higgs boson and one VFf* 1 , we get 



(A.154) 



A. 8 Higgs Bosons and Fermions 

For the cross-sections between one KK Higgs boson and one KK <Si7(2)w-singlet fermion, 
we have 

- /G) = 32 ^^ 2 (m 2 L + S /3) , (A.155) 
2 2 

crimtm^gt) = -J^L, (2L + 3/3) , (A.156) 

a(G+t m - tG+) = 2887Tm l s2p2c2 (l2c 2 w e 2 m 2 y 2 + m 2 L{+12c^(m 2 - s)y 2 e 2 

+ 9c 4 w (m 2 - s)y 2 t + 4se 4 } + s/3{4se 4 - 9c*mV}) , (A. 157) 

^(#i/ra - /G) = <7(Gfi/ fl i /#) , 
= v(Gtf R1 - /G+) , 

= <r(Gif m ^fG-), 
aiHitm 5 t) = <7(fTit L i fli) , (A.158) 
= a{Git R i -> flit) , 
= cr(G^t ra -> fl&) , 

= ^(G+6 L1 flt) . 

For the cross-sections between one KK Higgs boson and one KK 5C7 (2)vi/-doublet fermion, 
we get 

e 2 {Tfc w g 2 - 2 gi s w Yf) 2 , „ 
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a{ Hl ft - FG+) = 64n ^ 2s(32 (m 2 L + 8 0) , (A.160) 

- 1; wisy 2 ^ + «) + ( - 2l + °® ' (A - i6i) 

<x(G+t L1 ^+) = - (A.162) 



l l^-x e 2 (c w 52 - 2gis w Y f ) 2 , 2 , 
a^fcx-ftG )= 128W44s/32 + 

+ 32^ppjg2 (j/t^{*e 2 (4 - 2^) (A.163) 



-"W W 

+ c l s K s - m2 )Vt - s 2 w c 2 w sf3 2 y 2 t }) , 
where Tj? denotes the fermion isospin. 

a{H x b L1 - Gf) = \a{G+t L1 - , 

a(ffi/Li - G?/) = <7(Gi/ L i - #/) , 
= <r(Gff L1 -+ G ± f) , 
<7(fli/+ - /-G+) = a^i/f - /+G-) , 

= <7(Gi/+ - /"G+) , (A.164) 
= aCGfi/f - /+G~) , 
= a(G+/f - f+G) , 

where / stands for any lepton or quark, except tj_,\ and bn, and / + (/ _ ) denotes isospin 
+ 1/2 (isospin —1/2) fermions. 
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